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Abstract 

We propose a simple but useful parametrization of the flavor composition of 
ultrahigh-energy neutrino fluxes produced from distant astrophysical sources: 
(j) e : (p^ : <p T = sin 2 £ cos 2 ( : cos 2 £ cos 2 ( : sin 2 (. We show that it is possible 
to determine or constrain £ and £ by observing two independent neutrino flux 
ratios at the second-generation neutrino telescopes, provided three neutrino 
mixing angles and the Dirac CP-violating phase have been well measured in 
neutrino oscillations. Any deviation of £ from zero will signify the existence 
of cosmic v T and V T neutrinos at the source, and an accurate value of £ can be 
used to test both the conventional mechanism and the postulated scenarios 
for cosmic neutrino production. 
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I. INTRODUCTION 



High-energy neutrino telescopes are going to open a new window on the Universe, since 
they can be used to probe and characterize very distant astrophysical sources [1]. The 
promising IceCube neutrino telescope [2] , which has a kilometer-scale detector, is now under 
construction. If the relative fluxes of ultrahigh-energy v e (v e ), (v ) and v T (V T ) neutrinos 
are successfully measured at IceCube and other neutrino telescopes, it will be possible to di- 
agnose the relevant cosmic accelerators (e.g., their locations and characteristics) and examine 
the properties of neutrinos themselves (e.g., neutrino mixing and leptonic CP violation). 

Indeed, robust evidence for neutrino masses and lepton flavor mixing has been achieved 
from the recent solar [3], atmospheric [4], reactor [5] and accelerator [6] neutrino oscillation 
experiments. Due to neutrino oscillations, the neutrino fluxes observed at the detector 
$ D = {0^,0^,0^} are in general different from the source fluxes $ = {<f) e ,(j) fl ,(j) T }. Note 

that our notation is = + 6^ (for a = e,u,r), where 6^ and <b^ denote the 
z/ a -neutrino and I7 Q -antineutrino fluxes, respectively. The relation between <p u (or (f) v ) and 
4>u (or 4% p ) is given by 

< = E (K p J) > 

a 

^ = E (K P <#) ' 

a 

in which P a/3 and P a/3 stand respectively for the oscillation probabilities P{y a — > v£) and 
P(V a — > Vj). As the Galactic distances far exceed the observed neutrino oscillation lengths, 
P af3 and P af3 are actually averaged over many oscillations and take a very simple form: 

Pap = P aP = jl\Vj 2 \V pi \\ (2) 
i=\ 

where V ai and (for a, f3 — e, /j,,t and % = 1,2, 3) are just the elements of the 3x3 neutrino 
mixing matrix V. Eqs. (1) and (2) lead us to a straightforward relation between Q and 0°: 

$ = E {<t>aPap) ■ (3) 
a 

This relation indicates that the observation of <3> D at a neutrino telescope can at least help 1 

• to determine or constrain the flavor composition of cosmic neutrino fluxes [8,9], if three 
neutrino mixing angles and the Dirac CP-violating phase hidden in P aj3 have been 
measured to a good degree of accuracy (e.g., a precision of 10% or smaller relative 
error bars [10]); 



1 One may also use neutrino telescopes to test the stability of neutrinos [7,8], possible violation of 
CPT symmetry [8] and other exotic scenarios of particle physics and cosmology. 
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or 

• to determine or constrain one or two of three neutrino mixing angles and the Dirac 
CP- violating phase [11,12], provided the production mechanism of ultrahigh-energy 
neutrinos at an astrophysical source (e.g., the conventional source to be mentioned 
below) is really understood. 

Hence neutrino telescopes will serve as a very useful tool to probe both the high-energy 
astrophysical processes and the intrinsic properties of massive neutrinos. 

This paper aims at a determination of the flavor composition of cosmic neutrino fluxes 
at the source with the help of neutrino telescopes. Different from the previous works (see, 
e.g., Refs. [7-9,11,12]), our present study starts from a generic parametrization of the initial 
neutrino fluxes: 

{<f> e , </V , <t> T } = {sin 2 ^ cos 2 C, cos 2 f cos 2 C, sin 2 C} O , (4) 

where £ G [0, 7r/2] and ( G [0, 7r/2], and O denotes the total flux (i.e., the sum of three 
neutrino fluxes). Provided ultrahigh-energy neutrinos are produced by certain astrophysical 
sources (e.g., Active Galactic Nuclei or AGN) via the decay of pions created from pp and 
P7 collisions, for instance, their flavor content is expected to be 

K : K ■ tr} = {|:|:0} • (5) 

This "standard" neutrino flux ratio corresponds to £ = and tan £ = 1/ \[2 (or equivalently 
£ ~ 35.3°) in our parametrization. It turns out that any small departure of ( from zero will 
measure the existence of cosmic v T and V T neutrinos, which might come from the decays of 
D s and BB mesons produced at the source [13]. On the other hand, any small deviation of 
tan 2 £ from 1/2 will imply that the conventional mechanism for ultrahigh-energy neutrino 
production from the AGN has to be modified. Similar arguments can be put forward for 
the neutrino fluxes from some other astrophysical sources, such as the postulated neutron 
beam source [14] with 

{</> e : ^ : T } = {1 : : 0} (6) 
(or equivalently {£, (} = {vr/2, 0}) and the possible muon-damped source [15] with 

{</> e : ^ : T } = {0 : 1 : 0} (7) 

(or equivalently {£,C} = {0)0})- Therefore, we are well motivated to investigate how the 
true values of £ and ( for a specific astrophysical source can be determined or constrained by 
use of the second-generation neutrino telescopes and with the help of more precise data from 
the upcoming long-baseline neutrino oscillation experiments. This goal is indeed reachable, 
as we shall explicitly demonstrate in the remaining part of this paper. 

The remaining part of this paper is organized as follows. In section II, we derive the 
analytical relations between the neutrino flavor parameters (£ and () at an astrophysical 
source and the typical observables of neutrino fluxes at a terrestrial detector. Section III is 
devoted to a detailed numerical analysis of the dependence of those observables on £ and £. 
A brief summary of our main results is given in section IV. 
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II. OBSERVABLES 



Because of neutrino oscillations and the v T "regeneration" in the Earth [2], it is especially 
important to detect all three flavors of the cosmic neutrinos at a neutrino telescope. The 
sum of 0^, (f>® and (ff is equal to that of <f> e , <fi^ and <f> T , 



0o = 



+ </>„ + 



(8) 



as one may easily see from Eqs. (2) and (3). A measurement of O may involve large 
systematic uncertainties, but the latter can be largely cancelled out in the ratio of two 
neutrino fluxes. Therefore, let us follow Ref. [12] to define 



In b b\ = f g g & 



(9) 



as our working observables. We remark that these ratios are largely free from the systematic 
uncertainties associated with the measurements of 0j?, (j>® and 0°. In particular, it is 
relatively easy to extract R from the ratio of muon tracks to showers at IceCube [16], even 
if those electron and tau events may not well be disentangled. Since R e , R and R T satisfy 



R. 



+ 



R, 



+ 



Rr 



1 + iL l + R„ 1 + R 



1, 



(10) 



only two of them are independent. 

Note that Eqs. (6) and (7) represent two peculiar (non-standard) scenarios of cosmic 
neutrino production. In principle, one may also assume an exotic astrophysical source which 
only produces v T and V T neutrinos; i.e., {0 e : M : <p T } = {0 : : 1} or equivalently ( = n/2 
with unspecified £ in our parametrization. The expression of R a (for a — e, /x, r) can then 
be simplified in such special cases: 



Ra = < 



1-P„ 



P, 



l-P, 



P^ 



fin 



for {£,(} = {vr/2,0}, 
for {£,(} = {0,0}, 
for {£,C} = {*,tt/2}. 



(11) 



Even if the third case is completely unrealistic, it could serve as an example to illustrate the 
salient feature of R a defined above. 

Without loss of generality, we choose R e and R as two typical observables and derive 
their explicit relations with £ and (. By using Eqs. (4), (8) and (9), we obtain 



R„ 



P ee sin 2 £ + P^ cos 2 i + P Te tan 2 C 



sec 2 C - [P ee sin 2 £ + P^ cos 2 £ + P re tan 2 (] ' 
P eu sin 2 i + P uu cos 2 £ + P Til tan 2 C 



sec 2 C - [P e » sin 2 £ + P w cos 2 £ + P Tfl tan 2 (] 



(12) 
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The source flavor parameters £ and ( can then be figured out in terms of R e and R : 
r (p _p \_ r (p _p ) + (p P -P P 

2 e r Tfj, fifi I fx V re ^ie y 1 I MM re M e r M 



sin^ = 



( r _ p ) ( p _ p ") - (V - P ) (p - P ) 

V e re/ \ efi MM/ \ M T M/ \ ee M e / 



r fp _p")_ r fp _p") + fpp _pp") 

.2/- e \ MM e^y M V Me ee^ 1 ^ e/j ^e ee , , 

tan C = — j — r -j-^ 7 r— ^- , (13) 

where the notations r e = P e / (1 + P e ) and = P^/ (1 + P^) have been used to simplify the 
expressions. Indeed, r e = 0°/0 o an d r M = 0^/0o hold- One may in principle choose either 
(P e , R^) or (r e , r ) as a set of working observables to inversely determine £ and (. The first 
set has been chosen by a few authors (see, e.g., Refs. [10,12]) and will also be adopted in this 
paper. Note that the averaged neutrino oscillation probabilities P a/3 depend on three mixing 
angles (6* 12 , 9 23 , 9 13 ) and the Dirac CP- violating phase (5) in the standard parametrization of 
the 3x3 neutrino mixing matrix V [17]. Taking account of 9 13 < 10° [18] but 9 12 m 34° and 
9 23 ~ 45° [19], we express P af3 as the first-order expansion of the small parameter sin6> 13 : 

P ee = l-isin 2 2# 12 , 

P efl = - sin 2 29 12 cos 2 9 23 + - sin A9 l2 sin 2# 23 sin 9 13 cos 5 , 

P er = \- sin 2 29 12 sin 2 9 23 — j sin A9 12 sin 29 23 sin # 13 cos 5 , 

11 1 
Pnn = ^ ~ 2 S ^ n2 ^ 23 ~ 2 S ^ n2 ^ 12 C ° s4 ^ 23 ~~ 2 ^ 12 ^ 23 C ° s2 ^ 23 ^ 13 C ° S ^ ' 

11 1 

P = - sin 2 2# 2 3 sin 2 29 l2 sin 2 2 9 23 + - sin 4# 12 sin 49 23 sin 6 , 13 cos S , 

2 8 8 

P TT = 1 — - sin 2 26> 23 — - sin 2 2# 12 sin 4 #23 + - sin 46> 12 sin 26> 23 sin 2 #23 sin 9 13 cos 5 , (14) 

in which the higher-order terms, such as (9(sin 2 9 13 ) < 3%, have been safely neglected. One 
can see that the sensitivity of P af3 to 5 is suppressed due to the smallness of sin 9 13 . Hence the 
dependence of R a on 5 is expected to be insignificant. In addition, it is hard to distinguish 
between the cases of 9 13 = 0° and 5 = 90°, because it is the product of sin# 13 and cos 5 that 
appears in the analytical approximations of P a/3 . 

A global analysis of current neutrino oscillation data [19] yields 

30° < 9 12 < 38° , 
36° < 9 23 < 54° , 
0° < 9 l3 < 10° , (15) 

at the 99% confidence level. The best-fit values of three neutrino mixing angles are 9 12 ~ 34°, 
#23 ~ 45° and 9 l3 « 0° [19], but the CP-violating phase S is entirely unrestricted. Although 
the present experimental data remain unsatisfactory, they can be used to constrain the 
correlation between the parameters (£, C) an d the observables (R a) Rp). We shall illustrate 
our analytical results by taking a few typical numerical examples in the subsequent section. 
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III. ILLUSTRATION 



First of all, let us follow a rather conservative strategy to scan the reasonable ranges 
of (#12^23)^13)^) an d (£>C) so as to examine the sensitivities of (R e) R T ) to these six 
parameters. We take 5 G [0°, 180°] in addition to the generous intervals of three mixing 
angles given in Eq. (15), and allow £ to vary in the region £ G [0°, 90°]. As the amount of v T 
and V T neutrinos produced at those realistic astrophysical sources is expected to be very small 
or even vanishing, we restrict the parameter ( to a very narrow domain £6 [0°, 18°], which 
corresponds to sin 2 ( < 0.1 [13]. It should be noted that we do the numerical computation 
by using the exact expressions of P af3 instead of Eq. (14). Our results are shown in FIG. 1. 
Two comments are in order: 

1. The source parameter £ is quite sensitive to the values of the neutrino flux ratios R a 
(for a = e,fi,r). Even if three neutrino mixing angles involve a lot of uncertainties 
and the Dirac CP-violating phase is entirely unknown, a combined measurement of 
(R e , R^) or (R e , R T ) can constrain the value of £ to an acceptable degree of accuracy. 
This encouraging observation assures that the second-generation neutrino telescopes 
can really be used to probe the initial flavor composition of ultrahigh neutrino fluxes. 
Note that the standard pion-decay source {0 e : <f)^ : <f) T } = {1/3 : 2/3 : 0} will produce 
{0° : 0° : 0°} ~ {1/3 : 1/3 : 1/3} or equivalently R e ~ R^ ~ R T ~ 0.5 at the 
detector, as shown in FIG. 1, if 9 23 ~ 45° is fixed and 9 13 < 10° is taken. Such an 
expectation cannot be true, however, when £ deviates from its given value £ = 35.3° 
and (or) when 9 23 departs from its best-fit value 9 23 = 45°. More precise neutrino 
oscillation data will greatly help to narrow down the (R a ,£) parameter space. 

2. In contrast, the source parameter ( seems to be insensitive to R a (for a = e,fi,r). 
The reason for this insensitivity is two-fold: (a) the values of ( have been restricted 
to a very narrow range (0° < ( < 18°); and (b) the numerical uncertainties of three 
neutrino mixing angles and the Dirac CP- violating phase are too large. Provided 9 12 , 
9 23 , 9 13 and 5 are all measured to a high degree of accuracy in the near future, it will 
be possible to find out the definite dependence of R a on ( for a given value of £. 

To be more explicit, we are going to consider three typical scenarios of cosmic neutrino fluxes 
and illustrate the sensitivities of R e , R^ and R T to the source parameters (£, () and to the 
unknown neutrino mixing parameters (9 13 ,5). 

We argue that the simple flavor content of ultrahigh-energy neutrino fluxes from the 
standard pion-decay source could somehow be contaminated for certain reasons: e.g., a 
small amount of u e , and v T and their antiparticles might come from the decays of heavier 
hadrons produced by pp and jry collisions [13]. Similar arguments can also be made for 
the postulated neutron beam source and the possible muon-damped source, as our present 
knowledge about the mechanism of cosmic neutrino production remains very poor. Following 
a phenomenological approach, we slightly modify the scenarios listed in Eqs. (5), (6) and (7) 
by allowing the relevant £ and ( parameters to fluctuate around their given values. Namely, 
we consider 

• Scenario A: 30° < £ < 40° and 0° < ( < 18°, serving as a modified version of the 
standard pion-decay source (originally, £ = 35.3° and ( = 0°); 
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• Scenario B: 80° < £ < 90° and 0° < ( < 18°, serving as a modified version of the 
postulated neutron beam source (originally, £ = 90° and ( = 0°); 

• Scenario C: 0° < £ < f0° and 0° < ( < 18°, serving as a modified version of the 
possible muon-damped source (originally, £ = 0° and ( = 0°). 

For simplicity, we fix 9 12 = 34° and 9 23 = 45° in our numerical analysis. We take four typical 
inputs for the unknown parameters 9 l3 and 5: (a) 9 13 = 0° (in this case, 5 is not well-defined 
and has no physical significance); (b) 6 13 = 5° and 5 = 0°; (c) 6 13 = 5° and 5 = 90°; and (d) 
6 13 = 5° and 5 = 180°. Our numerical results for the sensitivities of R e , and R T to £ and 
( in scenarios A, B and C are shown in FIGs. 2, 3 and 4, respectively. Some discussions are 
in order. 

1. Scenario A in FIG. 2 . The neutrino flux ratios R e , R^ and R T are all sensitive to the 
small deviation of £ from its standard value £ = 35.3°. In contrast, the changes of 
three observables are very small when ( varies from 0° to 18°. This insensitivity is 
understandable: about half of the initial v., and V „ neutrinos oscillate into z/ and 
V T neutrinos, whose amount dominates over the survival amount of initial v T and V T 
neutrinos at the detector. In other words, 0^ ^> T and 0|? ~ 0° ~ 0° hold, implying 
that R a must be insensitive to small T or equivalently to small (. It is therefore 
difficult to pin down the value of ( from this kind of astrophysical sources. As pointed 
out in the last section, the result of R a in the 9 13 = 0° case (solid curves) is almost 
indistinguishable from that in the S = 90° case (dotted curves). The sensitivity of R a 
to 5 is insignificant but distinguishable, if the value of 9 13 is about 5° or larger. 

2. Scenario B in FIG. 3 . The neutrino flux ratio R e is sensitive to the small departures 
of £ and £ from their given values £ = 90° and ( = 0°. This salient feature can be 
understood as follows. Since v (or V ' ) and v T (or V T ) neutrinos at the detector mainly 
come from the initial v e (or V e ) neutrinos via the oscillation, the sum of 0^ and <ff is 
expected to be smaller than or comparable with the survival u e (or V e ) flux 0^. The 
roles of e and r are important in 0^ and 0° + 0°, respectively. It turns out that 
R e depends, in a relatively sensitive way, on £ through 0^ in its numerator and on ( 
through 0j^ + 0^ in its denominator. Note also that the nearly degenerate results of 
R e for 9 13 = 5° and 5 = (0°, 90°, 180°) are primarily attributed to the fact that P ee in 
the numerator of R e is actually independent of 5. On the other hand, the discrepancy 
between R e (9 13 = 0°) and R e (9 13 = 5°) in FIG. 3 results from the C(sin 2 6» 13 ) terms 
of P aj3 , which have been neglected in Eq. (14). In comparison with R e , the neutrino 
flux ratios R^ and R T are not so sensitive to the small changes of £ and (. But the 
measurement of R^ and R T is as important as that of R e , in order to determine the 
flavor composition of ultrahigh-energy neutrino fluxes at such an astrophysical source. 

3. Scenario C in FIG. 4 . In this case, in which 0° and 0j? mainly come from the initial 
M via v — > v e and — > v T oscillations, one can similarly understand the numerical 
behaviors of R e , R^ and R T changing with the small deviation of £ from its given value 
£ = 0°. To explain why three observables are almost independent of the fluctuation 
of (, we take 9 23 = 45° and 9 13 — > 0°, which guarantee P Ta = P (for a — e, /i, r) to 
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hold in the leading-order approximation. We can then simplify Eq. (12) in the £ — > 0° 
limit and arrive at the following result: 



oW " ^(-(p^+p^c) - i-p m ' (16) 



where ( is completely cancelled out. Therefore, the tiny dependence of R a on ( 
appearing in FIG. 4 is just a natural consequence of the 0(sin# 13 ) or £>(sin 2 # 13 ) 
corrections to Eq. (16). 

Although the numerical examples taken above can only serve for illustration, they do give us 
a ball-park feeling of the correlation between the source parameters (£, () and the neutrino 
mixing parameters (9 12 , 23 , 13 , 5) via the working observables (R e , R , R T ) at a neutrino 
telescope. This observation is certainly encouraging and remarkable. 



IV. SUMMARY 

We have proposed a simple parametrization of the initial flavor composition of ultrahigh- 
energy neutrino fluxes generated from very distant astrophysical sources: e : 0^ : T = 
sin 2 £ cos 2 C : cos 2 £ cos 2 ( : sin 2 (. The conventional mechanism and the postulated scenarios 
for cosmic neutrino production can all be reproduced by taking the special values of (£,()■ 
Of course, such a parametrization is by no means unique. An alternative, 

ill *Z> 1 X t , . 

^ ■ ^ ■ ^ = TTt : ITT : TVt (17) 

with x G [0, 1] and t G [0, oo) in general, is also simple and useful. For a realistic astrophysical 
source, it should be more reasonable to take < t <^ 1. Comparing between Eqs. (4) and 
(17), one can immediately arrive at x — sin 2 £ and t = tan 2 (. The (x, t) and (£, () languages 
are therefore equivalent to each other. 

After defining three neutrino flux ratios R a (for a — e, //, r) as our working observables 
at a neutrino telescope, we have shown that the source parameters £ and ( can in principle 
be determined by the measurement of two independent R a and with the help of accurate 
neutrino oscillation data. The standard pion-decay source, the postulated neutrino beam 
source and the possible muon-damped source have been slightly modified to illustrate the 
sensitivities of R a to small departures of £ and ( from their given values. We have also 
examined the dependence of R a upon the smallest neutrino mixing angle 9 13 and upon the 
Dirac CP-violating phase 5. Our numerical examples indicate that it is quite promising to 
determine or constrain the initial flavor content of ultrahigh-energy neutrino fluxes with the 
second-generation neutrino telescopes. 

How to measure R a to an acceptable degree of accuracy is certainly a big challenge to 
IceCube and other neutrino telescopes. A detailed analysis of the feasibility of our idea for 
a specific neutrino telescope is desirable, but it is beyond the scope of this paper. Here we 
remark that our present understanding of the production mechanism of cosmic neutrinos 
needs the observational and experimental support. We expect that neutrino telescopes may 
help us to attain this goal in the long run. 
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FIG. 1. Allowed regions of the neutrino flux ratios R a (a = e, fi, r) versus the source parameters 
£ and £, where we have scanned the 99% C.L. intervals of three neutrino mixing angles and taken 
the Dirac CP-violating phase 5 G [0°, 180°]. The horizontal and vertical lines in the (R a ,g) plots 
correspond to £ = 35.3° and R a = 0.5, respectively. 
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FIG. 2. Numerical illustration of scenario A , where # 12 = 34° and # 2 3 = 45° have typically 
been input in our calculations. The horizontal and vertical lines in the (£,R Q ) plots correspond to 
R a = 0.5 and £ = 35.3°, respectively. 
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FIG. 4. Numerical illustration of scenario C , where 9 12 = 34° and # 2 3 = 45° have typically- 
been input in our calculations. 
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